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Abstract 

We provide rigorous bounds for the error of the adiabatic approximation of quan- 
tum mechanics under four sources of experimental error: perturbations in the ini- 
tial condition, systematic time-dependent perturbations in the Hamiltonian, coupling 
to low-energy quantum systems, and decoherent time-dependent perturbations in the 
Hamiltonian. For decoherent perturbations, we find both upper and lower bounds on 
the evolution time to guarantee the adiabatic approximation performs within a pre- 
scribed tolerance. Our new results include explicit definitions of constants, and we 
apply them to the spin-1/2 particle in a rotating magnetic field, and to the supercon- 
ducting flux qubit. We compare the theoretical bounds on the superconducting flux 
qubit to simulation results. 

1 Introduction 

Adiabatic quantum computation [7] (AQC) is a model of quantum computation equivalent to 
the standard model [1] . A physical system is slowly evolved from the ground state of a simple 
system, to one whose ground state encodes the solution to the difficult problem. By a physical 
principle known as the adiabatic approximation, if the evolution is done sufficiently slowly, 
and the minimum energy gap separating the ground state from higher states is sufficiently 
large, then the final state of the system should be the state encoding the solution to some 
problem. 

The biggest hurdles facing many potential implementations of a quantum computer are 
the errors due to interaction of the qubits with the environment. However, the effects of such 
errors are different in AQC than in standard quantum computing. AQC is robust against 
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dephasing in the ground state, for instance [5], and some have suggested that noise in some 
regimes might actually assist adiabatic quantum computation [9]. 

The physical principle underlying AQC is established by the Adiabatic Theorem (AT). 
The AT bounds the run-time of the algorithm using the minimum energy gap of the system 
during the evolution. The AT itself has been recently subject to controversy [TH [231 121]) 
and cannot be applied directly to systems with noise or decoherence. There have been some 
numerical studies of AQC in the presence of noise OEj, and an analytic random- matrix study 
[T7] . Several recent studies have focused on the adiabatic approximation in open quantum 
systems using the density operator formalism [191 12H [23 El 122] • However, it is difficult to 
derive rigorous bounds with this approach because the dynamics involve a non-Hermitian 
operator without a complete set of orthonormal eigenstates. 

The subject of this work is the study of the adiabatic theorem in the presence of noise, 
perturbations, and decoherence. We consider several ways to extend the application of the 
theorem. Our statements, derived from Avron's version of the theorem [2], include explicit 
definitions of constants, so that we may apply them on examples. 

Experimental error for quantum computing experiments can be conveniently divided into 
three categories [2"5] : 

1. Coherent errors, due to a systematic implementation error such as miscalibration in a 
magnetic field generator. 

2. Incoherent errors, due to deterministic qubit-level differences in the evolution such as 
those caused by manufacturing defects. 

3. Decoherent errors, which are are random qubit-level errors due to coupling with the 
environment. 

In Section [2j we prove several extensions of the adiabatic theorem to handle these different 
types of error. For coherent errors, we provide a theorem for perturbations in the initial state 
of the system, and a theorem for systematic time-dependent perturbations in the Hamilto- 
nian. In the case of decoherent errors, we provide two new theorems, one for open quantum 
systems and one for noise modeled as a time-dependent perturbation in the Hamiltonian. 

In Section [31 we apply the new theorems to the spin-1/2 particle in a rotating magnetic 
field, a standard example for controversy regarding the Adiabatic Theorem [H [23], [23 [13] . We 
show our theorems make correct predictions about the error of the adiabatic approximation. 

Finally, in Section [4] we apply the new theorems to the superconducting flux qubit [15j . 
which has been proposed for adiabatic quantum computation [TT]. We use our theorems to 
determine a range of evolution times where the adiabatic approximation is guaranteed to 
perform well for a typical set of physical parameters and an apparently reasonable physical 
noise source. This provides the experimentalist with analytic tools for determining param- 
eters to guarantee the adiabatic approximation works well, without the need to perform 
numerical simulations. 
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2 Adiabatic Theorems for Noisy Hamiltonian Evolu- 
tions 



We begin with a Hamiltonian evolution Ti(s) parameterized by s £ [0, 1]. If we define r to 
be the total evolution time, then the Hamiltonian at time t is 7i(t/r). Thus, as r grows, 
TC(s) describes a slower evolution. Assume T~t(s) has countable eigenstates {IV'jts))} and 
eigenvalues A (s) < Ai(s)..., and consider the subspace 

*(a) = SpaQ{|^ Tn (s)),...,|Vn(s))} , (1) 

for some < m < n. Then the adiabatic approximation states that if the state of the system 
is contained in \l/(0) at t = 0, then at time t = s/r the state is contained in \I/(s). Notice 
that while the ground state \ipo(s)) may be important for physical reasons, the definition 
above allows consideration of a more general set of states. 

It is convenient to define an operator that computes the error of the adiabatic approxi- 
mation for a Hamiltonian evolution. We will need the projection operator P(s) that projects 
a state onto ^f(s). We will also need the unitary evolution operator U (s), that is the solution 
to Schrodinger's equation in the form 

U{s) = -itH(s)U(s) . (2) 

To compute the error of the adiabatic approximation, we apply P(0) to obtain the component 
of the initial state contained in \l/(0), evolve it forward in time by applying U(s), and then 
apply I — P(s) to compute the component of the state outside ^(s). For convenience, we 
define Q(s) = I — P(s), so the error operator is Q(s)U(s)P(0). 

In fact it will be most useful to bound the 2-norm of this operator, denoted \ \Q(s)U(s)P(0)\\. 
The 2-norm of an operator A is the square root of the largest eigenvalue of A^A, and in this 
case yields a bound on the magnitude of the output state, given a normalized input state. 

The version of the adiabatic theorem that we use to bootstrap our proof is based most 
closely on that of Reichardt [16] , which is based on that by Avron [2] (with later corrections 
[3] The differences between our theorem and Reichardt's theorem are 

• Our version of the theorem includes an explicit definition of constants, necessary to 
obtain quantitative bounds. 

• Our version of the theorem applies to subspaces rather than only a non-degenerate 
state. 

• We also present an integral formulation which provides better bounds when the energy 
gap is small for a very brief interval. 

Throughout the paper, we use units where h = 1. 

Theorem 2.1 (The Adiabatic Theorem (AT)). Assume for < s < 1 that H.{s) is twice 
differentiate, and let 

H(s) <h(s), H(s) <b 2 (s). (3) 
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Further assume that T~C(s) has a countable number of eigenstates, with eigenvalues 

^o( s ) < and that P(s) projects onto the eigenspace associated with the eigenvalues 

{A m (s),A m+ iO),...A n (s)}. Define 

,„M_\M \ ( a\ r x / min{A n+ i(s)-A n (s),A m (s)-A m _i(s)} m>0 
27/; C si 

d( s ) = i + — g( s ) = /-p(,). (4) 

Finally, assume 7(3) > a// s. T/ien we have 

.8£) 2 (0)6i(0) 8D 2 (s)6!(s) 



||g( s )c/( s )p(o)|| <- 



T7 2 (0) r7 2 (s) 

Jo r 7 2 (0 V 7lH / 
Proof. See Appendix A. □ 

Notice that the first two terms in Equation (jSJ) do not go to zero as s — > 0, which is a 
consequence of simplifications that were made to determine this bound. However, since AQC 
is the intended application of our results, we are only interested in the error bound at the end 
of the evolution, namely s — 1. Also, we will usually assume there are b\ > b\ (s), b 2 > 62(5), 
7 < 7( s ); an d > f° r s ^ [0, 1]- Then we can find a constant upper bound for the 

integrand in Equation (jSj) and thus bound the integral, resulting in the simpler expression 



T7 Z \ 7 



||Q( S )f/( S )P(0)|| < — - ( 2h + sb 2 + s ° K \ > 1 ) . (6) 



In fact, we will usually be interested in the AT for non-degenerate ground states, in which 
case m = n = and D = 1, and we can use the inequality 

8 / 1 (]b 2 \ 

\\Q(s)U(s)P(0)\\ < — 26x + sb 2 + s— 1 . (7) 
H \ 7 / 

Also notice our statement of the AT is consistent with the common interpretation of the 
theorem: if r ^> I/7 2 then the error in the adiabatic approximation is small. 



2.1 Coherent or Incoherent Errors 

Coherent or incoherent errors, due to systematic or deterministic perturbations, may occur in 
one of two ways: either as a perturbation in the initial condition or as a smooth perturbation 
in the Hamiltonian. In this section, we see how such errors affect the adiabatic approximation 
for a non-degenerate ground state. 

Let us first consider a perturbation in the initial state, 

|0(O))=r/(hMO)> + 5|<M) , (8) 
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where rj" 2 = 1 + \5\ 2 is a normalization factor, |^>o(0)) is the ground state of 7Y(0), and \4>±) 
is some state orthogonal to V'o(O)- It is not sufficient here to define the error of the adiabatic 
approximation as the norm of the operator Q(s)U(s)P(0), where P(s) is the projection onto 
\ipo{s)), since this does not depend on the initial state. The component of the final state 
which lies outside the ground state at normalized time s is Q(s)U(s)\<f>(0)), and so here we 
take this as the error. 

Theorem 2.2 (AT for Error in the Initial State (AT-Initial)). Let Ti(s) have the properties 
required by the AT, and let the initial state |0(O)) be as in Equation (TJj). Then the error is 
bounded as 

||Q( S )f/( S )|0(O)>|| < M (\S\ + (2h + sb 2 + s^-)) . (9) 



T7^ 
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Proof. Using the AT and the triangle inequality for operator norms, and noting that the 
norm of unitary and projection operators is unity, we have 



\\Q{8)u{ 8 )\m)\\ = \\Q(s)u{s) v {\Mo)) + s\<f>±))\\ 

= \\r 1 (Q(s)U(s)P(0)\MQ))+SQ(s)U(s)\ ( j )1 _) 
<\rj\(\\Q(s)U(s)P(0)\\ + \5\) 

166?' 



< \V\ \6\ 



T7 2 



2bi + sb 2 + s- 
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(10) 

(11) 
(12) 

(13) 
□ 



Now suppose there is a smooth perturbation in the Hamiltonian caused by a systematic 
error, so that 

H e (s) = H{s) + eA(s) . (14) 
Then we can use the AT on 7i e by observing that 



n e {s) 


< 


H{s) 


+ e 


A( S ) 


He{s) 


< 


H{s) 


+ e 


K{ S ) 



(15) 
(16) 



However, we must account for the difference in ground state between Ti e (s) and 7i(s). Since 
we want to measure error from the intended eigenstates of the system, not the perturbed 
eigenstates, the error operator is Q(s)U e (s)P(0), where we introduce the following notation: 

U e (s) The solution to U e (s) = —ir?i e (s)U e (s) . 

P e (s) The projection operator onto the ground state of H e (s). 

Q e (s) J-P e (s). 

7 e The minimum energy gap between the ground state and first excited 
state of Ti e (s). 
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Theorem 2.3 (AT for Systematic Error (AT-Error)). Assume that Ti e {s) has the properties 
required by the AT, and let 



n £ (s) 



<6i 



< bo 



(17) 



^/l-|(^o(O)|0o(O))| 2 = 5o 



)/l-|(Vo(l)l^(l)>| 2 = *i 



where |^o(s)) is the ground state ofH e (s) and \<f>o(s)) is the ground state ofH(s). Ifje > 0, 
then we have 



|Q(l)f/ e (l)P(0)|| < — 

T7 2 



Proof. We know from AT that 



26i + b 2 



166? 



+ S + 5i + 5 5i 



< 



26i + b 2 + 



but we want to find ||Q(1)[/ £ (1)P(0)||. So define 

AP(s) = P e (s) - P(s) . 

Then we have 

Q(l)Z7 e (l)P(0) = (Q e (l) + AP(1)) 17.(1) (P(0) - AP(0)) 

=g e (l)Z7 e (l)P e (0) - Q £ (l)f/ £ (1)AP(0) + AP(l)f/ £ (l)P £ (0) 
- AP(l)[/ e (l)AP(0) . 

Now, the 2-norm of unitary and projection operators is unity, so 

16fof 



HQ(l)^e(l)P 



< 



T7 £ 2 



26i + b 2 



le 



+ ||AP(0)|| + ||AP(1)|| + ||AP(0)||||AP(1) 
It remains to find ||AP(s)||. We hope to write 

\Ms)) = M(s)\Ms)) 



(19) 



(20) 

(21) 

(22) 
(23) 



(24) 



(25) 



for some unitary transformation M(s) that is close to the identity provided ipo{s) and <j>o{s) 
are close to each other. We use the Givens rotation, where the first basis state is |0o) an d 
the second is the complement of the projection of |^o) onto the first basis state: 



ei = \M S )) 



C2 



[i-\M*)){M*)\) \M*)) 



l - \{Ms)\M*))f 



(26) 
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The remaining basis states are chosen arbitrarily so long as the resulting basis is orthonormal 
and spans the Hilbert space. In that basis, Equation (1251) is realized by 



( C*{s) S{s) 
-S(s) C(s) 



S(s) 





( 1 \ 







V 1 J V J V° 

where we define 

c(s) = (Ms)\Ms)} 
s(s) = ^Ji-\(Ms)\Ms))\ 2 ■ 

We see that 

M(s)P e (s)M\s) = M(s)\Ms))(Ms)\M ] (s) 

= \Ms))(Ms)\ 

= P(s). 

Letting E(s) = I — M(s), we have 

AP(s) = P £ (s) - P(s) 

= M\s)P(s)M(s) - M\s)M(s)P(s) 
= M\s) [P(s),M(s)} 
= M\s)[E(s),P(s)\ . 

But we know E(s) and P{s): 

( 1-C*(s) -S(s) \ 
S(s) 1 - C(s) 



E(s) 



(\ 




P(s) 



so 



/ S(s) 
S(s) 



[E(s),P(s)} 



\ 



(27) 



(28) 
(29) 

(30) 
(31) 
(32) 

(33) 
(34) 
(35) 
(36) 



, (37) 



(38) 
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Finally, 



||AP0O|| = ||P e (*)-P(*)|| (39) 
= \\[E(s),P(s)}\\ (40) 
= S(s) . (41) 

Combining Equations (jUJ) and (T51]) yields the theorem. □ 

When it is inconvenient to compute 5q and 5i exactly, they can be bounded using the 
"sin(9) theorem" [2Ql page 251]: 



, < e||A(0)|| e||A(l)|| 

where Aq(s) is the energy of the ground state of TC e (s). If \q(s) is difficult to find, we can 
use the Bauer-Fike theorem [201 P a g e 192] to get 



* " 7(0) -t l|A(0)H S ^j(l)- t \'\A(l)r (43) 

where 7(5) is the energy gap of the unperturbed Hamiltonian H(s). 

A remarkable feature of AT-Error is that it does not depend directly on the magnitude 
of the perturbation term eA(s) except at the endpoints. It does not matter which path we 
take through state space, so long as we begin and end near the correct Hamiltonians and do 
not accumulate too much error along the way. 



2.2 Decoherent Errors 

Now we consider decoherent errors induced, perhaps, by noise in the environment. We first 
consider noise modeled as a coupled quantum system where the environment Hamiltonian is 
independent of time, and then as a classical time- dependent perturbation in the Hamiltonian. 

For the environment Hamiltonian 7i en v and interaction Hamiltonian eA(s), we can write 
the combined Hamiltonian / H e {s) as 

H £ (s) =H(s)®I + I®H env + eA(s) . (44) 

Direct application of the AT yields a very pessimistic result because the ground state of the 
composite system has, in the weak coupling limit, both the target system and the environ- 
ment in the ground state of their respective Hamiltonians. The target system remaining in 
the ground state and the environment tunneling to its first excited state will be considered 
a failure of the adiabatic approximation. An experimentalist probably cannot achieve the 
environmental ground state, and the energy gap between environment states is likely quite 
small so that the AT produces a very large error bound. 

One way to resolve this problem in the interpretation of the adiabatic approximation is 
to work with the density operator that results from the partial trace [19] : 

p(s) = Tr ew (p e (s)) , (45) 
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where p e (s) is the density operator associated with the state of the composite system 7i e {s). 
Usually we can write 

p(s) = C(s)p(s) , (46) 

where C(s) is a linear operator but not generally Hermitian. We might try to use this 
differential equation to prove an adiabatic theorem restated in terms of the expectation 
(0 o (s)|p(s)|0o( s ))> where |0o( s )) is the ground state of H(s). The problem is that C(s) does 
not have a complete set of orthonormal eigenstates, which is of great assistance in proving 
the AT. A rigorous bound on the error of the adiabatic approximation has yet to be found 
using this approach [T9l [281 l8l [2T1 [22] . 

The density operator approach sums together the set of states in the composite system 
whose measurement on the system of interest yields the ground state. Instead, below we will 
simply consider that set of states a subspace, and identify conditions where the usual AT for 
evolution of a subspace applies. 

For adiabatic quantum computation, we expect the energy gap 7 to be significantly larger 
than the temperature k^T . When the system is significantly coupled to only a small number 
N of nearby particles, then the range of relevant environmental energy levels is on the order 
of Nk B T, so 1 1 He™, 1 1 is also order NksT. If ||7i em ,|| < 7, we may use the AT. 

More generally, we must determine the error operator for evolution in the composite 
system. The projection operator of interest projects states in the composite system onto 
those states whose measurement reveals the original system to be in the ground state. If 
P(s) is the projection onto the ground state of H(s), then this operator is P(s) <S> I- Its 
complement is Q(s) ® / = / <S> I — P(s) <g> I, since Kronecker products distribute. Then the 
error operator is (Q(s) (g) /) U e (s) (-P(O) <g> I). 

Theorem 2.4 (AT for Coupling to Low- Temperature Environment (AT-Env)). Suppose we 
are given 

H € {s) =H{s)®I + I®H env + eA{s) , (47) 
and suppose we can choose w so that 

||Wen«||+2e||A(s)||<«;<7, (48) 

where 7 is the minimum energy gap between the ground state and first excited state oJT~L{s). 
Assume that 7i e (s) has the properties required by the AT, and assume that 7i env has M states 
and its ground state has zero energy. Let 



HJs) 



HJs] 



< h , (49) 



e 


|A(0)|| 


7 - 1 




— e | 


|A(0)|| 



d o = - — n^j — n 1 1 a / r\\ 1 1 > 01 = - — n^n — n m A/1M1 > ( 50 ) 



e 


A(l)|| 


7 - I 


| Heriv 


|| — e 


|A(1)|| 



•• < : 'i • Mr*:::!- (5i) 
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Then we have 

||(Q(1) ® /) 17.(1) (P(0) s 7)|| < ^ + 6 2 + 8(1 + D ^ ) +S + S 1 + ^ , (52) 

Tit \ It J 

•where r is the total evolution time. 

Proof. For e = 0, we can ignore 7i env and this theorem is simply the AT. So let us consider 
e > 0. We will do this by considering e > as a perturbation of the e = case. 

For e = 0, the eigenstates of TC £ (s) are simply the eigenstates of 77(s) tensored to the 
eigenstates of Ti, env , and the energy of those states is the sum of the energy of the state in 
77(s) and the energy of the state in H env . 

Define the ground state energy of 77(s) as A (s), the energy of the first excited state as 
Ai(s), and 7(5) = Ai(s) — A (s). Recall that the M energies of the TC env states are non- 
negative and less than 7. Then the first M eigenstates of H e (s) are the ground state of H(s) 
tensored with different eigenstates of 77 em) , and the rest of the states are some excited state 
of Ti(s) tensored with an environment state. 

In particular, the M th state of 77 e is the ground state of 7i tensored with the most 
energetic state of TC env , and thus has energy A (s) + ||77 em) ||. The M + 1 state is the first 
excited state of H tensored with the ground state of H env , and has energy Ai(s). So the 
energy gap between the first M states and the rest of the spectrum is is exactly ■y(s) — \ \H env \ \ . 

For positive e, these eigenstates are perturbed. Using the Bauer-Fike theorem [201 page 
192], we see that the gap is reduced by at most 2e ||A(s)|| in the presence of coupling, so the 
gap is still at least 7e- 

What we want is the adiabatic approximation of the evolution of the subspace formed 
by these M eigenstates, with a spectral width at most w and an energy gap of 7 e . Following 
the proof of AT-Error, define 

AP(s) = P e (s) -P(s) ®7 , (53) 

then we have 

(Q(l) <g> 7) U e (l) (P(0) ® 7) = (Q e (l) + AP(1)) U e (l) (i> (0) - AP(0)) (54) 

=g e (l)f/ e (l)P(0) - Q e (l)f/ e (l)AP(0) + AP(l)f/ e (l)P(0) 
- AP(l)t/ e (l)AP(0) , (55) 

We can bound 1 1 AP(s) 1 1 using the fact that the singular values of AP(s) are given by the 
sines of the canonical angles between P e (s) and P(s) <8>7 [201 page 43], the "sin(O) theorem" 
[20| page 251], and the Bauer-Fike theorem [20], page 192]: 



e 


|A( 


»l 




7(s) - 1 


H 


3nn | 




C 


|A(a)|| 



AP(s) < — '' , /M — - , (56) 



so 

AP(0) < 5 AP(1) < 5 1 . (57) 
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Now we are ready to apply the AT: 



||(Q(i)® /) u £ (i) (P(o)®/)|| < ||g e (i)t/ e (i)p e 

8P 2 



+ S + 5\ + S 5i 



(5c 



< 



T7 £ 2 



26i 



8(1 + D)6? 



7c 



+ 5 + 5i + 5 5 1 . (59) 



□ 



Now let us consider another model of decoherent noise, namely a time-dependent per- 
turbation in the Hamiltonian. There is a problem applying the AT directly, because the 
time-dependent perturbation is a function of true time t, not the unitless evolution param- 
eter s. So as r grows, more noise fluctuations are packed into the interval s G [0, 1], causing 
||dW/ds|| to diverge. Then there is no bound b\ greater than ||d7i/ds|| that is independent 
of r. In fact this problem was the source of confusion in the recent controversy surrounding 
the adiabatic theorem [T3 | 1231 123] . 

We will need to consider Hamiltonians T~C T {s) that depend on both s and t. We define 
the following notation: 

U T (s) The solution to U T (s) = —iTH T (s)U T (s) for a fixed r. 
-P-r(s) The projection operator onto the ground state of Ti T {s). 

Q T (S) I-P T (S). 

7 T (s) The energy difference between the ground state and first excited state of 7t T (s) 

Theorem 2.5 (Adiabatic Theorem for Hamiltonian Evolutions on Two Time Scales (AT-2)). 
Suppose, for any fixed r, that 7i T (s) has the properties required by the AT. Further assume 
there are real functions gi(r) and g-i (r) such that 



H T (S) 

for all t. If there is a 7 m j 



<9i(r) 
so that 



Hr{s) <<? 2 (r) 



< 7„ 



< 



for all s and t, then we have 

\\Q T (s)U T (s)P T (0)\\ < 



2gi (t) + sg 2 (r) + s 



7mm 



(60) 



(61) 



(62) 



Proof. The lemma we are trying to prove is the union of special cases of the AT, when the 
AT is applied to one-parameter projections of the original Hamiltonian. 

For fixed r, we consider TC T (s) as a one-parameter Hamiltonian to which the usual AT 
will apply. Then by the AT, we can write 



||Q T (s)tf T (s)P T (0)|| < 



2<7i(r) + sg 2 (r) + 



In 



But we can do this for any r, so the lemma holds. 



(63) 
□ 
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Now we can apply AT-2 to the case where there is an evolution performed on some scaled 
time s, with an additive noise Hamiltonian 7i no i se (t) that is a function of real time t = st: 



H T (s) = H(s) + H noise (sT) . 



(64) 



We define the error operator for the noisy Hamiltonian as Q(s)U T (s)P(0). The projection 
operators refer to the unperturbed Hamiltonian because success should be defined in terms 
of the intended states. 

Theorem 2.6 (Adiabatic Theorem for Noisy Hamiltonian Evolutions (AT-Noise)). Suppose 
for any fixed r, that H T (s) = 7~C(s) +Ti. no i se (sT) has the properties required by the AT. Assume 



< Ci 

< di 



dt 2 



d 2 



< C 2 

< d 2 



yJi-\(MQ)\M0))\ 2 = s 



y/i-\(Mi)\Mi))\ 2 = 5i 



(65) 
(66) 
(67) 



where \ipo(s)) is the ground state of 7i T (s) and |0o(s)) is the ground state ofTi(s). Further 
assume that there is a ^j no i se so that 



< In 



(68) 



for all s and r . Then we have 
||Q(l)C/ r (l)P(0)|| < 



16d?\ , / 16ci 
d 2 + — -± r + 2d! 1 + - 



Tnoise L V 

+ 5 + Si + 5 5 1 . 



In 



2ci + c 2 



16c? 



Proof. Evidently 



ds ds dt 



d 
ds 



H T (s] 



< ci + rd\ 



d 2 



< c 2 + r 2 d 2 ■ 



T 

(69) 



(70) 
(71) 

(72) 
(73) 



Substitution of Equation ( JTTj) and Equation ( 1731 into AT-2 yields, for all s and r, 

16c? 



\\Q T (s)U T (s)P T (0)\\ < 



7 2 ■ 

/ noise 



16d? \ , / 16ci 

d 2 + — -± r + 2d! 1 ' 



In 



In 



2ci + c 2 



In 



(74) 
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As in the proof of AT-Error, we define 



AP(0) = P T (0) - P(0) AP(1) = P T (1) - P(l) . (75) 

Then for s = 1 we have 

Q(1)U T (1)P(0) = (Q T (1) + AP(1)) C/ r (l) (P T (0) - AP(0)) (76) 
=Q T (l)f/ r (l)P r (0) - Q T (l)C/ r (l)AP(0) + AP(l)[/ r (l)P T (0) 

- AP(l)£/ r (l)AP(0) . (77) 

Now we bound the norm of the error: 

16d? \ _ , ( . 16ci \ / n 16c? \ 1 

rf 2 + - — L r + 2di 1 + - — - + 2ci + c 2 + 



l)C/ r (l)P(0)|| < 

i noise L \ ] noise/ \ I noise / \ In 

+ ||AP(0)|| + ||AP(1)|| + ||AP(0)|| ||AP(1)|| . (78) 
Using the Givens rotation just as in the proof of AT-Error, we have 

||AP(0)||=5 ||AP(1)||=5 1 , (79) 

which, when substituted into Equation (1781) . completes the proof. □ 
Several observations can be made about this result: 

1. As with AT-Error, if it is inconvenient to compute Sq and 5\ exactly, they can be 
bounded using the "sin(O) theorem" combined with the Bauer-Fike theorem [2"Ul page 
192]: 

r ^ | |7"^n.otse(0) 1 1 r ^ I \7~t-noiseiX) 1 1 (Rfi) 

- 7 (o) - ||7Wo)ll 1 " 7(1) - \\n noise (i)\\ ' 1 ] 

where j(s) is the energy gap between the ground state and first excited state of H(s). 
Also, S and Si can be taken as zero if H-noisei®) = W no i se (r) = 0. In general, we expect 
them to be quite small if H no i se (t) is several orders of magnitude smaller than TC(s). 

2. When r is small, the 1/r term dominates. This term is exactly the bound from the 
(noiseless) AT. It shows there is always a positive lower bound on the running time of 
the adiabatic algorithm to guarantee a particular error tolerance. 

3. When r is large, the first term dominates. In fact, we can see that in the presence of 
noise, there is always some sufficiently large r beyond which the adiabatic approxima- 
tion may perform poorly. So given an error tolerance, there is always an upper bound 
on the running time for the adiabatic algorithm, beyond which the theorem cannot 
guarantee the tolerance to be met. 

4. If there is a great deal of noise, and thus di is large, the constant term (with respect to 
r) could become as large as 0(1) and there could be no running time for the adiabatic 
algorithm which results in an accurate calculation. 
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We are also interested in a lower bound on the error of the adiabatic approximation in 
the presence of noise. A lower bound could be used to prove that a certain amount of noise 
was unacceptable for AQC, because it would guarantee failure of the algorithm for some 
level of noise. However, it will be difficult to get a non-trivial lower bound, since there are 
time-dependent perturbations which yield zero error in the adiabatic approximation, better 
than might exist without the perturbation. To see this, define 

H A (s)=H(s)+i/T[P(s),P(s)}, (81) 

where the term i/r[P(s), P(s)} is the perturbation. Avron proved [2J, as do we in the 
appendix (Theorem lA.lj) . that the evolution of TIa(s) satisfies the intertwining property, 
assuming Ti(s) is non-degenerate, countably dimensional, and twice differentiable. The 
intertwining property can be written 

U A (s)P(0) = P(s)U A (s) , (82) 

where U A (s) is the unitary operator associated with TCa(s), and P(s) is the projection onto 
the ground state of 7i(s). Then 

Q(s)U A (s)P(0) = Q(s)U A (s)P(0) (83) 
= Q(s)P(s)U A (s) (84) 
= , (85) 

since Q(s)P(s) = 0, so the adiabatic approximation is perfect for H(s) if the perturbation 
term i/r[P(s), P(s)] is added to it. Notice further that the perturbation gets arbitrarily 
small as r grows. 

Finally, we also observe that noise that commutes with the Hamiltonian does not cause 
any state transitions, because it has no effect on the eigenstates - in other words, it causes 
no coupling between states. For instance, consider the Hamiltonian on N particles 

N 

H(s) = M(s)J2<rh ( 86 ) 
i=i 

where M(s) is a real scalar function representing a time-dependent applied magnetic field. 
Noise in the magnetic field M(s) results in a perturbation that commutes with TC(s), and 
has no effect on the error of the adiabatic approximation. 

3 Application to the Spin- 1/2 Particle in a Rotating 
Magnetic Field 

Recently Tong et al. [22] presented an example of a Hamiltonian evolution for which the 
adiabatic approximation performs poorly. The Hamiltonian is for a spin-1/2 particle in a ro- 
tating magnetic field. Here we apply AT-Noise to their example. Their evolving Hamiltonian 
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IS 



7~t(t) = — ^ (a x sin 6 cos ut + o y sin 6 sin ut + cr z cos #) (87) 
which we represent in the z basis as 

^ (t) ""TVe^sinfl -cos^ ) " (88) 

Suppose # is small. We can think of the time-independent diagonal component of the 
Hamiltonian as the intended Hamiltonian, and the wobbling off-diagonal component as a 
noise term operating on an independent timescale. 

The eigenstates of 'Hit) depend on t, but the eigenvalues do not. So the energy gap is 
constant and in fact equal to u . Thus one might think that the adiabatic approximation 
works well, predicting that a particle starting out in the spin-down state stays in the spin- 
down state under this Hamiltonian evolution. 

We will see below that if the wobble is at a resonant frequency with respect to the 
energy difference between the spin-up and spin-down states, the wobble induces a complete 
transition from the spin-down to spin-up state. So the adiabatic approximation eventually 
fails in the most complete sense possible in this example. However, we will also see that the 
AT-Noise correctly provides an increasing error bound with time, because the s-derivatives 
in this example increase with r. 

We can rewrite the Hamiltonian using t = st as 

_ujo[ cos 9 e-^ ST sm9\ 
nr[S> 2 \e iuJST sm9 -cos 9 J ' 1 } 

We can also compute the first two derivatives: 

d horsing / -ie"*"" \ 

HrS = S -iusT n 7 I 90 ) 



ds w 2 V ie 

and 



d 2 „ , . , u 2 u q t 2 sin 9 ( 
dT 2HAs) = 2 { e 

We can compute the norms of these matrices exactly, giving 



Hl(s)H T (s) = ^ 



4 



1 
1 



dH T (s) 
ds 



dH T (s) 
ds 



u 2 u 2 T 2 sm 2 ( 1 



1 



(92) 
(93) 
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and 



Thus we can write 



d 2 H T (s) 
ds 2 



d 2 n 7 



ds 2 



< ^ (1 + \u sin 6\ t + uj 2 |sin 6>| r 2 ) 



cu^t 4 sin 2 9 ( 1 





1 



(94) 



(95) 



for all < s < 1 and for k — 0,1, 2. Also, 7 r (s) = uj q for any s, r [23], so we let ^ no i se = ujo- 
Schrodinger's equation can also be solved exactly for the Hamiltonian in Equation ( 1891) . 
Define u) = \/ uj 2 + uj 2 + 2ujquj cos 6. From Tong et al. [23] the unitary time evolution opera- 
tor for this system is 



U T {t) 



(COS (f ) + z ^ocose gin (ffltj j 

in (— ) 



(f)e- 



■iu;t/2 



(cos (f ) - i^fsal sin (|*)) e **t/2 



Therefore the error operator for the adiabatic approximation is 

Q(t)U T (t)P(0) = 



)U (t) f 1 ° 
1 i rl j I 











i^af^sin (f ) e iuJt / 2 



so 



a> sin 6* 




sin 




to 





\\Q(t)U T (t)P(0)\\ 

If the perturbation is resonant, then uj = —uj cos8 so uj = \ujq\ sin#. Then we have 

' uj Sm^t" 



(96) 

(97) 
(98) 

(99) 



\Q(t)U T (t)P 



sin 



(100) 



As an example, assume that 9 = 0.001, uj = 10, c^o = —10. Let xi T ) be the error bound 
defined by the adiabatic theorem. Then we can calculate ||-P r (s) — P(s)|| exactly to get 
5 = 5i = 0.0005, and so we have 



and 



X ( T ) = 0.00900025 + 0.04r 



\\Q(s)U T (s)P(0)\\ = |sin(0.005sr) 



(101) 



(102) 



Figure [T] illustrates our result. Not only is the bound consistent with the true error, it 
has the same qualitative behavior, increasing linearly with r. 
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Figure 1: A plot of the error bound x( T ) f° r parameters 9 = 0.001, u> = 10, and u = —10, 
compared to the true error for the Tong et al. [23] example. We can see that the adiabatic 
approximation gets worse as r gets larger, as observed by Tong et al. [23J. However, our 
error bound from AT-Noise remains valid. 

4 Application to a Superconducting Flux Qubit 

Next we apply AT-Noise to the superconducting flux qubit of Orlando et al. [15], proposed 
for use in adiabatic quantum computation [11] . With this qubit, the adiabatic evolution may 
be as simple as monotonically varying an applied magnetic field. 

Consider the four-junction qubit shown in Figure[2] We will follow the analysis of Orlando 
et al. [J5]. The dynamical variables are the phases 0j across the four Josephson junctions, 
however flux quantization in each loop gives us two constraints: <pi — <p2 + 4>3 = ~ 27r/i and 
04 — 03 = — 2vr/ 2 , where f\ and f'2 are the magnetic frustrations in each loop. So there are 
two degrees of freedom, which we define as <p p = (0i + 02 ) /2 and <fr m = (0i — 02 ) /2, where 
fa = f2, and fb = fi + ^2/2. Then the Hamiltonian can be written as 

n —&M-^W m + u{M > (103) 

where M p and M m are constants, and the potential U(<f) p , m ) is defined as 

U((f> p , <p m ) = Ej[2 + 2{3-2 cos(0 p ) cos(0 m ) - 2(3 cos(vr/ a ) cos(27r/ b + 20J] , (104) 

where Ej and j3 are constants. 

At fi = f2 = 1/3, we have fb = 1/2 and f a = 1/3, and U(4> p ,(p m ) has minima, or wells, 
at 4> p = and <p m = ± cos~ 1 (/3/2), symmetric about the P axis. By varying fi and f'2, we 
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Figure 2: (Adapted from Orlando et al. [15]) The circuit schematic for the superconducting 
flux qubit. The X's represent Josephson junctions, and the main qubit loop is formed by 
junctions Jl, J2, and J3. Junction J4 allows tuning of the effective properties of junction J2 
through control of the frustration f 2 - 

can tilt the potential so that one well is deeper than the other, and we can adjust the barrier 
height. We can approximate the Hamiltonian with a two-state system. A Hamiltonian 
evolution that begins at the degeneracy point and varies f\ can can be written: 



where t\ and r\ are parameters that can be estimated with the WKB approximation. For 
the qubit parameters recommended by Orlando et al. [15], r\ = A.8Ej and t\ = 10~ 3 Ej, 
where Ej is a constant. A typical value for Ej is 200 GHz • h = 1.256 THz • h. We choose 
e = —.0002, so that the Hamiltonian changes from proportional to a x at s = to equally- 
weighted a x and a z terms at s — 1, because this seems a natural milestone in the evolution 
to the cr 2 -dominated final Hamiltonian. 

There are a couple of sources of noise in this qubit. One source of noise in a superconduct- 
ing flux qubit is noise in the critical current of the Josephson junctions [26], which decreases 
as 1/T where T is temperature. Such noise would result in variations in the weights of the 
terms in Equation f)104p . Another source is noise in the magnetic flux bias generated by 
nearby current- carrying wires on the chip. Current carrying wires could be used for nearby 
measurement devices or to perform a gate operation on the qubit with an RF pulse. Since we 
have from above a two-state Hamiltonian parametrized by flux bias, we consider this latter 
noise. Orlando et al. [15] estimated that a nearby wire of typical dimensions and carrying 
lOOnA current would cause a difference in either f\ or f 2 of Af = 10~ 7 . Let us assume that 




(105) 
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there is there is approximately 0.5 nA of noise on the wire introduced by the current source. 
Further suppose the power of the noise scales as inverse frequency \ jv from v min = 2.5 GHz 
to Vmax = 3.5 GHz, so that we include the qubit frequencies throughout the evolution. 

There are clever means of simulating \ jv noise with discrete models, such as by summing 
independent bistable ffuctuators (a.k.a. Random Telegraph Noise) [6]. However, this results 
in non-differentiable Hamiltonians so is not appropriate for us. Instead, suppose we want to 
write down a formula for a noise source with 1/u power spectrum in the range v m in to v max . 
Let n be an integer, with n = 100 in the following example. Define Az/ = {y max — v m in)/ n 
and Vj = v m in + J 'Ai/ for j = l...n. Then we can define two independent noise functions, 
representing variation in the magnetic frustration in our qubit, as 



(106) 
(107) 



where £ij and £ 2 ,j are phase factors chosen uniformly at random and C = 10~ 10 MHz^ 1 / 2 , 
chosen to agree with the 0.5 nA noise. The Hamiltonian for noise in the magnetic frustration 
is 



7-lnoise{t) = Ni{t)ri<r z + N 2 {t)(r 2 (T z - wa x 



108) 



where w = IAEj for the chosen qubit parameters. 

Evaluating the functions numerically over an interval much larger than the longest wave- 
length reveals the bounds \M(t)\ < 4.9100 • 10~ 10 , Mi(t) < 9.1100 • 10" 6 MHz, and 

Niit) < 0.1667 MHz 2 . 

Recalling that t = st, where s is unitless, we are ready to compute derivatives and norms 
of the whole Hamiltonian: 



H T (s) = H(s) + Hnoiseisr) 

= -t\(T x + sr 1 ea z + Af 1 (t)r 1 a z + J\f 2 (t) (r 2 a z - wa x ) 



H r (s) 



Ht{s) 



riea z + T^Ni{t)ria z + r^Af 2 (t)(r 2 a z - wa x ) 

< 1206.4 MHz + r • 84.7149 MHz 2 , 

d 2 d 2 
r2 -^W) ri a z + r 2 —N 2 (t)(r 2 o z - wa x 

< r 2 1.5502 ■ 10 6 MHz 3 . 



(109) 
(110) 

(111) 

(112) 

(113) 

(114) 



Observe that since s is unitless, T~C(s) and its s-derivatives all have units of energy, but since 
h — 1, energy units are inverse time units. 
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Figure 3: A plot of the error bound x( T ) f° r the superconducting flux qubit example. The 
crosses represent results of numerical simulations of the error. Since the noise commutes 
with the dominant term in the final Hamiltonian, the bound is a substantial overestimate. 
Nonetheless the qualitative shape between the bound curve and the simulation data agree. 

We also need to compute the minimum energy gap. In this case, it occurs at s — 0, and 
the energy gap is % ise = 2ti = 2513 MHz. 

Finally, we need to find 5q and 8±. We compute the projection operators directly and 
obtain the bounds Sq = 1.800 • 10 -6 and 5\ = 9.117 • 10~ 7 . From AT-Noise, we have 



where if r is supplied in microseconds, the resulting bound is unitless. 

This generates a hyperbolic curve with a vertical asymptote at r = and a linear 
asymptote for large r, shown in Figure [3) Recall this curve represents the norm of the error 
operator and its square represents the probability of error in this system. 

To check our results, we would like to compute the error of the adiabatic approximation 
numerically. However efficient numerical simulation of this system requires some care. A 
straightforward solution to Schrodinger's equation 



in the a z basis will have rapidly oscillating phases that make the solutions very unstable and 
time-consuming to compute. Instead, we will rewrite Schrodinger's equation for this system 
in a basis whose phase rotates with time. 





i-\m)=n T {t/T)\m) 
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To begin, we choose a time-dependent eigenbasis of Ti T {t/r) with the property 
{4> n (t)\4> n (t)) = 0. In other words, the Berry phase, also known as geometric phase, is zero. 

Lemma 4.1. There is a time dependent eigenbasis {\4> n {t))} with the property that 
(<t>n(t)\Mt)) = 0for aUnJfl. 

Proof. Suppose {\4> n {t))} does not have this property. Then for each n, define the variable 
£ n {t) as follows: 

£ n (t) =i [ (Mw)\Mw))dw (117) 



and let 

\a n (t)) = e*»®\<f> n (t)) . (H8) 
Then {\a n (t))} has the desired property. To see this, we simply compute it: 

(a n (t)\d n (t)) = (Mt)\e~ lUt) [iUt)e iUt) \4>n(t)) + e*»W|^(f)>] (119) 

= i£ n (t) + (Mt)\Mt)) (120) 

= i ■ l{<P n (t)\^ n (t)) + (0 n (t)|0„(t)) (121) 

= . (122) 

□ 

Let us write the solution \ip(t)) in terms of the basis states \(j> n (t)) with energies E n {t) as 
follows: 

m)) = E c -W e " J/ ° t£n( ^i^w) • ( 123 ) 

n 

Then in the case of two states, assuming the eigenstates are labeled in increasing order with 
respect to their eigenvalues, the norm of the adiabatic error operator is simply |ci(£)|. 

Let us substitute the representation in Equation (11231) into Schrodinger's equation (Equa- 
tion (I116P ). and left multiply by {<t> m (t)\. After simplification, this yields 

emit) = -^c„(t)e-^o(^W-^W)^ (0m(t) |^ (t)) (124) 

Evidently the n = m term in the sum is zero in the chosen basis. 
Since we have 

n T {t/r) = {-h -M 2 (t)s 2 )a x + + + -A/" 2 (*)r 2 J a z , (125) 

if we define 

a(t) = -h -N 2 (t)s 2 , (126) 
b{t) = t -^ll + M 1 {t)r l +N 2 {t)r 2 , (127) 

T 

WHcar'gg), (128 ) 
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then we can diagonalize the Hamiltonian easily in terms of a(t) and b(t). We choose the 
cotangent for numerical stability because 6(0) ~ 0. 

Eo{t) = -y/a?(t)+P(t) E 1 (t) = y/a*(t) + V*(t) (129) 

cos I -tt I / y sm 

Now, we would like to compute (<f)o(t)\<j)i(t)) and (0i(t)|0o(^))- 

_ s i n ( ffl] \ / _ I sin ( m.) 0( t ) 

<*,(*)i*.w> = ( ^/J-( ; cos( ^ (() i (i3i) 

(132) 

(0 1 (t)|0 o (t))= I . >'< I • I "! COS lfx! W I (133) 




cos 




sin | 


v 2 


6{t) 




2 





-isin(^)^) 



It remains to compute 9(t), which can be done with implicit differentiation. 



a 2 (t) a(t)b(t) - a(t)b(t) 
a 2 (t) + b 2 {t) tf(t) 
a(t)b{t) - a(t)b(t) 
a 2 {t) + b 2 {t) 



Finally, the equations of motion are 

c 



(134) 



cotW)) = ^ (135) 
a{t) 

-esc (9{t))6(t) = ^r-r (136) 

9(t) = sin 2 (6>(t)) A(t)6(t) -° (t) ^ (t) (137) 

a 2 (t) 



(138) 
(139) 



(t) = -^(i)e-»J?VWW*»® (140) 
c x (t) = ^(^e^/oV^H+^H^M (Ml) 



We provide these equations to a differential equation solver, ode23, in Matlab. Care must 
be taken with the integral in the exponent. We need not recompute the integral entirely at 
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each time; rather we cache the intermediate values of this integral. Thus at each evaluation 
we only integrate on the interval from the last cached time to the current time. 

This method was used to produce the numeric results in Figure [3J The parameters of the 
system are those previously described in the example. There are 571 noiseless data points 
and 125 noisy data points. A new set of random phases was generated for each noisy point, 
and each point took up to 5.5 CPU hours to compute. The workstation used had dual Xeon 
3.06 GHZ processors with hyperthreading enabled (thus four effective CPUs) and 6GB of 
RAM, running Red Hat Enterprise 3. 

The simulation data in Figure [3] is several orders of magnitude less than the bound. The 
fact that the bound is an overestimate is not surprising since the noise term commutes with 
the dominant term in the final Hamiltonian of the evolution. However the qualitative shape 
between the bound curve and the simulation data is the same, and the bound does provide 
a simulation-free guarantee of error for an interval of r. 

5 Conclusion 

We provide rigorous bounds for the adiabatic approximation under for four sources of ex- 
perimental error: perturbations in the initial condition, systematic time-dependent per- 
turbations in the Hamiltonian, coupling to low-energy quantum systems, and decoherent 
time-dependent perturbations in the Hamiltonian. 

We applied the new results to the spin- 1/2 particle in a rotating magnetic field, which is a 
standard example for discussing controversy in the adiabatic theorem [U [23], [27] . We showed 
that our theorem makes correct predictions about the error of the adiabatic approximation 
function of time. 

We also applied the new results to the superconducting flux qubit proposed by Orlando 
et al. [15], with time-dependent perturbations in the applied magnetic field. This qubit has 
properties that make it a candidate for quantum adiabatic computation [IT]. Because our 
version of the adiabatic theorem does not have unspecified constants, we are able to make 
numerical predictions about this qubit. We showed that for a particular amount of noise 
on superconducting wires near a qubit with ideal physical parameters, we could guarantee a 
small error in the adiabatic approximation provided that the evolution time was set within 
a particular interval. 
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Appendix A: Proof of the Adiabatic Theorem 



Our proof of the adiabatic theorem follows closely those by Avron et al. [2] (later corrections 
exist [3] [12]), Reichardt [16], and Jansen et al. [10]. The purpose of revisiting the proof is 
to have explicit definitions of constants. 

After reviewing some properties of projection operators, we will introduce the version of 
Schrodinger's equation that will be used and the assumptions that it requires. Then we will 
introduce some essential lemmas, and finally the proof. 

A.l Properties of projection operators 

Before embarking on the proof, it will be helpful to review some properties of projection 
operators. First, we will enumerate some elementary properties. Then, we will introduce 
the resolvent formalism for rewriting projection operators as a contour integral. 

Define the commutator [A, B] as [A, B] = AB — BA, and A = dA/ds. Let Ti.(s) be some 
Hamiltonian with countable eigenstates {(^(s))} and eigenvalues {Aj(s) : j > 0}. Let P(s) 
be the orthogonal projection operator onto the subspace 

*( S ) = Span{|^ m ( S )>,...,|^ n ,( S )>} , (A-l) 

for some < m < n. Thus 

k 

P(s) = Y;Ms))tA(s)\. (A-2) 

Let Q(s) = I — P(s) be the orthogonal complement of P(s). Then the following properties 
hold: 

1. P(s) =P 2 (s). 

2. P(s) = P(s)P(s) + P(s)P(s), obtained by differentiating both sides of Property ([I]). 

3. P(s)P(s)P(s) = 0, obtained by multiplying Property ([2]) from the left by P. 

4. Q(s)P(s) = 0, using the definition of Q and Property ([T|). 

5. P^(s) = P(s) and Q^(s) = Q(s), where > indicates the conjugate transpose. This is 
evident from Equation (1A-2|) . 

6. ||P(s)|| = ||Q(s)|| = 1. Recall that the norm of an operator P(s) is defined to be the 
maximum of ||P(s)|x)|| for choices of normalized states \x). For a projection operator, 
the maximal choice is a vector in the plane of projection, and in that case P{s) \x) = \x). 
So ||P(s)|| = 1. For Q(s), choose \x) orthogonal to the plane of projection. 
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7. [H(s), P(s)] = 0. To prove this, take some state \<f>) and rewrite it as |0) = 2^j>o Cj\ipj(s)). 
Then 

n(s)P(s)\ ( j ) ) = Y,Cjn(s)P(sM(s)) (A-3) 

i>o 

n 

= ^>W(a)|^00) (A-4) 

n 

= J2 c i X MMs)) (A-5) 



and 



P( a )W( a )|^> = £p( a )W(a)c,|^(a)>+ £ P(a)W(a)c,|Vi(a)> (A-6) 

j=m jt^[m,n] 
n 

= X)c,-P(«)A i («)|^(a)>+ £ c,P( s )A,( S )|^( S )) (A-7) 

j=m jt^[m,n] 
n 

= X)c J -A i ( a )|^(a)> . (A-8) 



j=m 



We will make use of the resolvent formalism to bound the projection operators. Define 
the resolvent of a Hamiltonian 7~£(s) to be 

P(z; W(s)) = (W(s) - ziy 1 . (A-9) 

Suppose we can draw a contour r(s) in the complex plane whose enclosed region includes 
the eigenvalues corresponding to ^(s) and excludes the rest of the spectrum of 7i(s). Then 
we can rewrite the projection operator P(s) in terms of a line integral of the resolvent 
R(s, z) = (T~l(s) — zl)~ l around this contour: 

P(s) = — — <p R(s, z)dz . (A-10) 
2tt« Jt(s) 



A. 2 Schrodinger's equation 



We rewrite Schrodinger's equation in terms of unitary evolution operators and the scaled time 
s, rather than using state vectors and real time t. Doing so will introduce the assumption 
that Ti(s) has a continuous, bounded second derivative. 

The usual expression of the time-dependent Schrodinger's equation is 

K^^=nt)\m)- (A-ii) 
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Setting t = st, \<f>(s)) = |^(sr)) = \i/)(t)), and TC T (s) = TC(sr) = H(t), we can substitute 
and apply the chain rule for derivatives to get 



h\j>(s)) = -irH T {s)\(t>{s)) , (A-12) 

where the dot indicates the s-derivative. It will be convenient in this proof to choose our 
units so that h — 1. Also we will assume that all subsequent state vectors, Hamiltonians, 
and time evolution operators are functions of the normalized time parameter s, so we can 
drop the subscript r from 7i T . Thus we will write 

\<j>(s)) = -iTH(s)\<j>(s)) . (A-13) 

Now define U(s) so that for any [(f)(0)), we have U(s)\(j)(0)) = \4>(s)) where \4>{s)) is the 
solution to this equation. Then we proceed as in [18]. Assume that H(s) has a continuous 
bounded derivative; then \4>(s)) has a continuous bounded second derivative. Thus the 
remainder for the first-order Taylor expansion is well-defined. For some point s* G [s, s+As], 
we get 

U(s + As)|0(O)> = \(f>(s + As)) (A-14) 

= |0( S )) + |0(s))As + |0( s *))^ (A-15) 

= |0(s)> - irH(s)\(f)(s))As + 0{As 2 ) (A-16) 
= U(s)\(/)(0)) - irn(s)U(s)\(f)(0))As + 0(As 2 ) . (A-17) 

Since this is true for any 1 0(0) ) we can write 

lim U{S + A A S) ~ U{S) = -rrH{s)U{s) (A-18) 

or, equivalently, 

U{s) = -itH(s)U(s) . (A-19) 



Equation flA-19j) is the form of Schrodinger's equation that we will rely on for the rest of the 



proof of the adiabatic theorem. 
A. 3 Essential Lemmas 

Recall that the adiabatic approximation states that a system of Hamiltonian T~C(s), initially 
in some state in \l/(0), evolves to approximately some state in \P(s) at time t = sr. To 
compute bounds on the error of this approximation, we will identify a Hamiltonian 71a that 
has exactly this property. Define 

H A (s) = H{s) + - \P{s), P(s)] (A-20) 

T L 



26 



where P(s) is the projection operator onto \&(s). Evidently Ha is a 1/r perturbation of H, 
where r is the scale factor between normalized time and unnormalized time. Define Ua(s) 
to be the unitary evolution operator that is the solution to Schrodinger's equation for Ha, 
namely 

U A (s) = -itHa(s)U a (s) . (A-21) 

The important property of Ha{s) can be restated as follows. If a system is initialized in 
^(0) at time s — 0, the state at time s under evolution by the Hamiltonian Ha{s) is entirely 
contained in We can write this property, known as the intertwining property, using 

P(s) and Ua{s) as defined in the previous paragraph. 

Theorem A.l (The Intertwining Property). For allO < s < 1, letH(s) be Hermitian, twice 
differentiate, non- degenerate, and have a countable number of eigenstates. Let Ua{s) and 
P(s) be defined as previously. Then 



U A (s)P(0) = P(s)U A (s) 



(A-22) 



Proof. Noticing that Ua(s) is unitary, we can rewrite the claim as P (0) = U a (s)P(s)Ua(s) 
Since Ua{0) = I this is certainly true for s = 0. So it is sufficient to show that 



d_ 

ds 



U\(s)P{s)U A {s) 



. 



Applying the product rule for derivatives we get 



d 
ds 



U\{s)P{s)U A {s) 



d_ 
ds 



U\(s)P(s) 



U A (s)P(s) 



U A {s) + U A (s)P(s)U A (s) 
U\(s)P( S ))Ua(s) + U\(s)P(s)U a (s) 



(A-23) 

(A-24) 
(A-25) 



Now observe that U A 
its matrix entries. Further, recall that Ua(s) = —itHa{s)Ua{s). So 

U\{s) = (Z7 A ( a ))t 

= (-itH a (s)U a (s)) ] 
= +itU\{s)H\{s) 
= itU\{s)H a {s) , 
since Ha is Hermitian. Substituting, we get 
d 



{Ua) since the derivative of a matrix operator is the derivative of 



(A-26) 
(A-27) 
(A-28) 
(A-29) 



ds 



U A (s)P(s)U A (s) 



U A (s)P(s) +irU A {s)H A {s)P{s)) U A (s) 
+ U a (s)P(s)(-z)tHa(s)U a (s) 
U A (s) (P(s) + ztH a (s)P(s) - trP(s)H A (s) ) U A (s) 



-U a (s)[P(s)+zt[Ha(s),P(s)})U a (s) 



(A-30) 
(A-31) 

(A-32) 
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Now we will work on the inner term. We use the properties that [H(s), P(s)] = 0, P(s)P(s)P(s) 
0, P 2 (s) = P(s), and P(s) = P(s)P(s) + P(s)P(s). 

[H A (s), P(s)] =H A (s)P(s) - P(s)H A (s) (A-33) 

= (h{s) + i (P(s)P(s) - P(s)P(s)fj P(s) 

- P(s) (h(s) + I (p( a )P(a) - P(a)P(a))) (A-34) 
= [H(s),P(s)] + -P(s)P 2 ( S ) - -P(s)P(s)P(s) 

T T 



- -P{s)P{s)P{s) + -P 2 (s)P(s) (A-35) 
r r 

:-P(s)P 2 (s) + -P 2 (s)P(s) (A-36) 
r r 

:-P(s)P(s) + -P{s)P{s) (A-37) 
r r 

:-P(s) . (A-38) 
r 



We can substitute this into the original expression to get 



d_ 

ds 



U\(s)P(s)U A (s) 



U A (s) (P(s) + ir[H A (s), P(s)]) ^(s) (A-39) 
171(a) (^(«) - ^(«)) Ua(s) (A-40) 



= . (A-41) 

□ 

Notice that this implies an intertwining property for the orthogonal complement: 

U A (s)Q(0) = Q(s)U A (s) . (A-42) 

In the proof of the adiabatic theorem we will make use of the twiddle operation. For a 
fixed s, let P be a projection operator onto and assume the eigenvalues corresponding 

to ^(s) are separated by a gap from the rest of the eigenvalues. Define 

X(s) = —l R(s,z)X(s)R(s,z)dz (A-43) 
2tt« Jr(s) 

where T(s) is a contour in the complex plane around the eigenvalues associated with the 
eigenstates onto which P(s) projects, whose enclosed region excludes any other eigenvalues 
of TC(s). We will need the following property of the twiddle operation. 

Lemma A. 2 (The Twiddle Lemma). Assume % = 1. For a fixed s, let P be a projection 
operator onto ^{s), and assume the eigenvalues corresponding to \&(s) are separated by a gap 
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from the rest of the eigenvalues. Define Q = I — P , and let X be a bounded linear operator. 
Then 

qxp = -q ( \n A ,x] -- 



p 



(A-44) 



Proof. We begin by observing that since P 2 = P and QP = 0, 

-Q[X, P]P = -Q(XP - PX)P 
= -QXP . 



(A-45) 
(A-46) 



Further, since the identity operator commutes with everything, [zl, R(z)XR(z)} = 0. Then 



2m 



[H,X] = Hi — : & R{z)XR{z)dz - — - i R{z)XR{z)dz H 



2m 



— I {HR{z)XR{z) - R{z)XR{z)H) dz 
2m J r 

— l\H,R{z)XR(z)]dz 
m J v 

[H - zI,R{z)XR{z)]dz 



2m 
1 

2m 



(A-47) 
(A-48) 
(A-49) 
(A-50) 



Now we use the fact that (TC-zI)R = I, that X does not depend on z, and Equation ( 1A-10I) 
to write 



[H, X] = — / (XR(z) - R(z)X) dz 
2m J r 



= -[X,P\. 

Also, using the definition of Ha in ( 1A-20I) . we have 

QXP = Q[X, P]P 
= -Q[H,X]P 

= -Q 



Ha, A 

All we need to finish the proof is to show 



r ll 



P,P 



X 



P 



(A-51) 

(A-52) 
(A-53) 

(A-54) 
(A-55) 

(A-56) 



Q 



P,P 



x 



P = Q 



P,X 



p . 



(A-57) 
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Using PPP = 0, P 2 = P, and P = PP + PP, we have 

Q[[P,P}X}P = Q[PP-PP,X]P (A-58) 

= Q (pP - PP) XP - QX (pP - Pp) P (A-59) 

= (I - P) (PP - PP) XP-(I - P)X (PP - PP) P (A-60) 

= PPXP - PPXP + PPXP - XPP + PXPP (A-61) 

= PXP - PPXP - XPP + PXPP (A-62) 

= (J - P)PXP - (I - P)XPP (A-63) 

= QPXP - QXPP (A-64) 

= Q[P,X]P. (A-65) 

□ 



A. 4 Completion of the Proof 



Now we are ready to prove the adiabatic theorem. The operator Q(s)U(s)P(0) takes the 
^(0) component of a state vector, and evolves it for (normalized) time s under T~C(s), then 
computes the error from the adiabatic approximation. The theorem guarantees an upper 
bound on the norm of this operator. 

Theorem A. 3 (The Adiabatic Theorem). Assume for < s < 1 that Ti{s) is twice differ- 
entiable, and let 



H(s) <6i(s) 



H(s) <h(s), 



(A-66) 



Further assume that Ti(s) has a countable number of eigenstates, with eigenvalues 

\q(s) < Ai(s)..., and that P(s) projects onto the eigenspace associated with the eigenvalues 

{X m (s),X m+1 (s), ...X n (s)}. Define 



w(s) = \ n (s) - \ m (s) , 7(s) 
2w(s) 



min{A n+ i(s) - A n (s), A m (s) - A m _i(s)} m > 

m = ' 



P> s = 1 



A n+ i(s) - A n (s) 

g( s )=/-p(s). 



7T7(s) ' 

Finally, assume 7(5) > all s. Then we have 

.8P 2 (0)6i(0) 8P 2 (s)&!(s) 



(A-67) 



||0(a)l7( a )P(0)|| <- 



T7 2 (0) r7 2 (s) 

s 8P 2 (r) /8(1 + P(r))6 2 (r) 



r~i 2 (r) 



7(r) 



+ &2(^) dr . 



(A-68) 



30 



Proof. By multiplying by the identity and applying Theorem IA. 1 1 (the intertwining property) , 
we can write: 



Q{s)U(s)P(0) = Q(s)U A {s)tf A (s)U(s)P(0) 
= U A (s)Q(0)U A (s)U(s)P(0) . 



(A-69) 
(A-70) 



Since ||£7a(s) 



= 1, if Q(0)U A (s)U(s)P(0) 
adiabatic approximation is small. In fact, if we define 

W(s) = U A (s)U(s 



is small the magnitude of the error in the 



(A-71) 



then W(s) satisfies a useful integral equation, and we prove the adiabatic theorem by bound- 
ing ||Q(0)lU(s)P(0)|| instead of working directly on Q(s)U(s)P(0). To find the integral 
equation, we need to compute W(s). Using the product rule for derivatives, Schrodinger's 
equation, and the definition of H A (s) in (1A-20I) . we have 



W(s) 



U A (s)U(s) + U A (s)U(s) 
-irU A (s)H(s)U(s) + irU A (s)H A (s)U(s) 



-U A (s) 



-U A (s) 



P(s),P(s) 
P(s),P(s) 



U(s) 

U A (s)W(s) . 



(A-72) 
(A-73) 
(A-74) 

(A-75) 



Clearly W(0) = I and so 

W(s) 



U A (r)[P(r),P(r)]U A (r)W(r)dr 



(A-76) 



It will be useful sometimes to refer to the kernel of this integral equation, so we define 



K(r) 



U\{r) 



P(r),P(r) 



U A {r) 



(A-77) 



Now we can use Equation flA-76j) to rewrite ||Q(0)W / (s)P(0)||. Using the fact that 
Q(0)P(0) = 0, we can write 



Q(0)W(s)P(0) 



Q(0)K(r)W(r)P(0)dr . 



(A-78) 



Our plan is to rewrite the integrand to obtain an expression where all but one term has 
a 1/r factor. Integration by parts on the remaining term will ensure all terms have a 1/r 
factor. Then we can factor out the 1/r and bound the operators in each term to yield the 
Adiabatic Theorem. 

To obtain this expression, we need to introduce a P(0) in the middle of Equation ( 1A-78I) 
so that we can apply Lemma IA.21 To do so, we will use the fact that Q(0) = Q{0) 2 to 
introduce another Q(0), and then use the fact that Q(0)K(r) = K(r)P(0). 
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To show that Q(0)K(r) = K(r)P(0), we use intertwining properties, the fact that 
Q(r)P(r) = 0, and the properties P 2 (r) = P(r) and P(r)P(r)P(r) = 0: 



Q(0)U A (r)[P(r),P(r)]U A (r) 


(A-79) 


U A (r)Q(r)[P(r),P(r)]U A (r) 


( A-l 




U A (r) 


( K Q{r)P{r)P{r) - Q(r)P(r)P(r)^j U A {r) 


(A-l 


11) 


U A (r)Q(r)P(r)P(r)U A (r) 


(A-i 


12) 


U\(r) 


(>(r)P 2 (r) 


- P(r)P(r)P(r)\ U A (r) 


(A-l 


13) 


U A (r) 


P(r),P(r) 


P(r)U A (r) 


(A-l 


14) 


U A (r) 


P(r),P(r) 


U A (r)P(0) 


(A-l 


15) 


K(r)P(0) . 




(A-l 


16) 



Then we can rewrite 



Q(0)W(s)P(0) 



Q(0)K(r)W(r)P(0)dr 
Q(0) 2 K(r)W(r)P(0)dr 
Q(0)K(r)P(0)W(r)P(0)dr . 



(A-87) 
(A-88) 
(A-89) 



Now we use the definition of K(r), the properties P 2 (r) = P(r) and Q 2 (r) = Q(r), and the 
intertwining property again: 

Q(0)W(s)P(0) = - [ Q(0) 2 U A (r) P(r),P(r) U A (r)P(0) 2 W(r)P(0)dr (A-90) 

Q(0)U A (r)Q(r) \P(r), P(r)l P(r)U A (r)P(0)W(r)P(0)dr . (A-91) 

We would like to apply Lemma [A. 21 with 



X(r) = P(r),P(r) 



(A-92) 



In order to apply the lemma, we need to show X(r) is a bounded linear operator. Clearly 
X(r) is linear, and since P(r) has unit norm then it is sufficient to show that P(r) has a 
bound. 

To bound the norm of P(r) we will use the resolvent formalism. We first need to bound 
the norm of the resolvent R(r, z). Notice that if the eigenstates of TC(r) are {\4>j(r)) : j > 0}, 
then 



R(r, z)=J2 



1 



i>o 



Xj(r) -z 



\^(r))(^(r)\, 



(A-93) 
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Figure 4: Visualization of the resolvent contour T(s). The eigenvalues of H(s) (represented 
with black dots) all lie along the real axis since T~i(s) is Hermitian. Notice T(s) lies at least 
r y{s)/2 from any eigenvalues. The length of T(s) is 717(5) + 2w(s) = D(s)7rj(s). Observe 
that D(s) is the ratio of the length of T(s) to the circumference of a circle of radius r y(s)/2. 



so the norm of R(r, z) equals the inverse of the minimum distance of z to an eigenvalue of 
T~L{r). So we need to choose the contour T(r) around the eigenvalues of \&(r) to maximize 
the minimum distance of T(r) to any eigenvalue. Also, to obtain the best bound on the path 
integral, we will want to minimize the length of T(r), given that maximum minimum distance. 
We choose T(r) consisting of two semicircles connected by lines, forming a pill-shape. The 
semicircles are centered at A m (r) and A n (r), and they have radius r y(r)/2. Figure 0] illustrates 
this choice, which bounds the norm of R(r, z) at 2/7(7") and the length at D{r)'K n f{r). 

We can check the tightness of this choice by using it to bound the norm of P(r), which 
we know is unity. We have 



\P(s) 



2m J r(r) 
= D(r) , 



(A-94) 

(A-95) 
(A-96) 

so the approximation is tight for D(r) = 1. When D(r) > 1, it is complicated by the fact 
that the closest eigenvalue is not always the same at different points on T(r). 

The elements of R(r, z) are rational functions of the elements of T~C{r), which are assumed 
to be differentiate. So we can apply the quotient rule for derivatives to determine that 
R(r, z) is differentiable for z not an eigenvalue of TC(r). 

We proceed by differentiating both sides of the equation 



R{r,z){H{r) - zl) = I 
and multiplying both sides by R(r, z) on the right. We thus obtain 

R(r,z) = -R(r,z)H{r)R(r,z) . 



(A-97) 



(A-98) 
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So by Equation (lA-lOj) 



P(r) = ^—.1 R(r, z)H(r)R(r, z)dz 
2m Jr(r) 



(A-99) 



Also, recall that 



Ti{r) < bi(r), so we can bound the norm of the integral in Equa- 
tion flA-99j) with a rectangle approximation. Using our formula for the length of T(r), we 
get 



~ 2tt K J n J 7 (r) 2 
2D{r)h{r) 



P(r) < ^D(r)wir)—^- (A-100) 

(A-101) 



j(r) 

Finally, we can bound ||X(r)||. Using the definition of X we have 



\X(r)\\= [P(r),P{r)\ (A-102) 

(A-103) 
(A-104) 



< 2 


P(r) 


\\P(r 


AD{r)b x 






7 (r) 





Thus we can apply Lemma [A. 21 We remove the extra Q(r) and P(r) the same way they 
were introduced, and use Schrodinger's equation. 



Q(0)W(s)P(0) = / Q(0)U A (r 



i 

T 



P(r),X(r) 



H A {r),X(r 
Q(0)U A (r)H A (r)X(r)U A (r)P(0)W(r)P(0)dr 
Q(0)U A (r)X(r)H A (r)U A (r)P(0)W(r)P(0)dr 
Q(0)U A (r) \P{r),X{r) \ U A (r)P(0)W(r)P(0)dr 
Q(0)U A (r)H A (r)X(r)U A (r)P(0)W(r)P(0)dr 
Q(0)U A (r)X(r)U A (r)P(0)W(r)P(0)dr 



U A (r)P(0)W(r)P(0)dr 



(A-105) 



Q(0)U{(r) 



P(r),X(r) 



U A (r)P(0)W(r)P(0)dr 



(A-106) 



Evidently the last two integrals have a 1/r factor, and we only need to work on the first 
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integral. We will integrate it by parts, using 

A(r) = X(r)U A (r)P(0)W(r) (A-107) 

dA = X(r)U A (r)P(0)W(r)dr + (x(r)U A (r) + X(r)U A (r)) P(0)W(r)dr (A-108) 

B{r) = U A (r) (A-109) 

dB = iTU ] A {r)H A {r)dr . (A-110) 
Applying the integration by parts to J dBA yields 



U A {r)H A (r)X(r)U A {r)P{0)W{r)dr = --U A (r)X{r)U A (r)P(0)W(r 

+ ~ [ U A (r)X(r)U A (r)P(0)W(r)dr 



r=0 



I 

T 

+ Z I U A (r)X(r)U A (r)P(0)W(r)dr 

s 



r 



s 



+ - I UUr)X(r)U A (r)P(0)W(r)dr . (A-lll) 
T Jo 

When we substitute, we see that the last integral cancels with the second integral in Equa- 
tion (IA-1061) . so we obtain 



Q(0)W(s)P(0) = --Q(0)U A (r)X(r)U A (r)P(0)W(r)P(0) 



r=0 



-Q(0) / U A (r)X(r)U A (r)P(0)W(r)P(0)dr 
■'" Jo 

+ l -Q(Q) J U\(r) (x(r) - [P(r),X(r)]) U A (r)P(0)W(r)P(0)dr . 

(A-112) 

To finish the proof, we need to bound each of the three terms on the right. We will do 
this by applying the triangle inequality to all the operators in Equation (1A-112|) . Unitary 
operators and projection operators have unit norm, and we have bounded P(r) already, so 

it remains to bound the norms of X(r), X(r), and W(r). As dependencies we will also need 
to find the norms of P(r) and X(r). 



P(r) 



To bound P(r) we need to compute R(r,z). Using the product rule for 



derivatives, 



-R(r, z) =4-R(r, z)H(r)R(r, z) (A-113) 
dr 

= (R(r, z)H(r)^j R(r, z) + R(r, z)H(r)R(r, z) (A-114) 

= (k{r, z)H{r) + R(r, z)H(r)) R(r, z) + R(r, z)H{r)R(r, z) (A-115) 

=R(r, z)H{r)R(r, z) + R(r, z)H(r)R(r, z) + R(r, z)H(r)R(r, z) . (A-116) 
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Since 



R(r) < \\R(r)\\ 2 H(r) < 46i(r)/7(r) 2 , we have 



R(r, z) 



< 



IQbAr) 2 Ab 2 (r) 



7(r) 3 7(r) 2 
A fAb^ 2 



7(r) 2 \ 7(7") 



So, following the reasoning used to bound 



P(r) 



P(r) 



2ni 



R(r, z)dz 



r(r) 



1 4 fAbAr) 2 



2D(r) /4&i(r 



7 ( r ) ^ 7 ( r ) 



+ &2(r) 



X(r) : By Equation (IA-430 . we have 



<^ / \\R(r,z)\\ \\X(r)\\ \\R(r,z)\\dz 

^ JT(r) 

1 , . . . 2 AD(r)bAr) 2 
- 2vr ,v ; v ; 7 (r) 7 (r) 7 (r) 
8D(r) 2 6i(r) 
7(r) 2 



X 



: Notice X(r) = P(r),P(r 



so 



x(r; 



[P(r),P(r)] 

<2 P(r) ||P(r 
4£>(r) /46i(r) 2 



7 ( r ) ^ 7 ( r ) 



+ & 2 (r) 



36 



4. 



X{r) 



We have 



X(r) 



1 

27 
1 

27 
1 

"27 
1 

"27 



— * i?(r, z)X(r)R(r, z)dz 
dr J T{r) 

I ^-R(r,z)X(r)R(r,z)dz 
Jtm dr 



(A-128) 
(A-129) 



r(r) 



R(r, z)X(r)R(r, z) + [R(r, z)X{r) + R(r, z)X{r) ) R(r, z)dz 



■R(r, z)X(r)R(r, z)H(r)R(r, z) - R(r, z)H(r)R(r, z)X(r)R(r, z) 

(A-130) 

Now since \\R(r, z)\ \ < 2/7 (r), we get 



r(r) 

+R(r,z)X(r))R(r, z)dz 



1 /I661W 4 

<-^(r)D(r) (^||X( r )|| + _ 



X(r) 



8D(r) 2 f8h(r 



7(r) 2 \ 7(r) 



(A-131) 
(A-132) 



5. 



^(r) : Recall from Equation (1A-75I) that 



W{r)=-U\(r) P(r),P(r) U(r)W{r) 



(A-133) 



We know that ||W(r)|| = U\{r) = \\U(r)\\ = 1, and remember that X(r) 



P(r),P{r) 



So we can apply the triangle inequality to get 



W(r) <\\X(r)\\ 

4D(r)6i(r) 



< 



The resulting bounds are 
2D(r)b 1 (r) 



P(r) 
W(r) 



< 



< 



7(r) 
AD(r)h(r) 
7 (r) 



X(r) 
X(r) 



< 



< 



7(r) 



8D(r) 2 6!(r) 

7(r) 2 
8D{r) 2 [Sb^rf 



(A-134) 
(A-135) 



7(r) 2 \ j(r) 



+ h(r] 



(A-136) 
(A-137) 
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Now let us apply these bounds to Equation (1 A- 112ft by taking the norm of both sides: 



\\Q(0)W(s)P(0)\ \ < - Q(0)U A (r)X(r)U A (r)P(0)W(r)P(0) 



r=0 



+ 



+ 



Q(0) / U\{r)X{r)U A {r)P(0)W{r)P(0)dr 
Jo 

Q(0) / U\(r)(X(r) + [P(r),X(r)])U A (r)P(0)W(r)P(0)dr 
Jo 



(A-138) 



We can further simplify this by noting that the norm of each integral is less than the 
integral of the norm of its integrand. Further, we use the triangle inequality and the fact 
that the norm of unitary operators and projection operators are unity: 



1 



\\Q(0)W(a)P(0)\\ <- X(0) 



+ 



X(s) 



X(r) W{r) 



< 8^(0)6i(0) 8D 2 (s)6i(s 



X(r) 



T7 2 (s) 



+ 



r 7 2 (0) 

8D 2 (r) /8(1 + D(r))b\(r) 



o ^7 2 (r) 

we get 



7(r) 



Finally, from Equation (j7 

||g( a )^(*)P(o)||<||g(o)w( a )p( 

8D 2 (0)6!(0) 8D 2 (s)6!(s) 



< 



T7 2 (s) 



[P{r),X(r)} dr) (A-139) 



bo(r) I dr . 



+ 



r 7 2 (0) 

8£> 2 (r) /8(1 + D(r))b\{r) 



to r-f 2 (r) V 7(0 
We also know that ||Q(s)£/(s)P(0)|| < 1 by the triangle inequality 



+ b 2 (r) I dr . 



(A-140) 
(A-141) 

(A-142) 
□ 
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